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We investigate effects of interactions between substitutional defects on the properties of defect 
crystals at low temperatures, where defect motion is governed by quantum effects. Both, thermal 
and dynamical properties are considered. The influence of interactions on defect motion is described 
via a collective effect. Our treatment is semiclassical in the sense that we analyze collective effects in 
a classical setting, and analyze the influence on quantized defect motion only thereafter. Our theory 
describes a crossover to glassy behavior at sufficiently high defect concentration. Our approach is 
meant to be general. For the sake of deflniteness, we evaluate most of our results with parameters 
appropriate for Li-doped KCl crystals. 



I. INTRODUCTION 

Interacting quantum impurities occur in various sys- 
tems, such as spin glasses, magnetic systems with random 
couplings, and ionic crystals with substitutional defects. 
Spin glasses have been studied as a prominent example 
for systems far from thermal equilibrium. Though real 
spin glasses are of quantum nature, most theoretical ap- 
proaches rely on classical models. Work on true quantum 
spin glasses has been mainly directed towards elucidat- 
ing the nature of the phase transition to a frozen low- 
temperature state and at studying the influence of quan- 
tum fluctuations on the order parameter (for a recent re- 
view, see Ref. 1). Likewise with quantum ferromagnets.^ 
Theoretical work has concentrated on the quantum phase 
transition that occurs at T = 0. In disordered quantum 
ferromagnets, additional interest stems from the fact that 
rare fluctuations of the disorder configurations leading 
to Griffiths singularities have a much stronger effect in 
the presence of quantum fluctuations than in the cor- 
responding classical systems.^. Pertinent experiments 
have been performed on the diluted dipolar Ising mag- 
net LiHocYi_cF4,'^ which exhibits both, magnetic and 
spin glass order, depending on the concentration c of Ho 
ions. Agreement between experiment and theory for this 
system is, however, still not in a satisfactory state. 

Early work on ionic crystals with substitutional im- 
purities, such as lithium or cyanide in various alkali 
halides (KCLLi, KBr:CN,...), focused on the thermal 
properties;^'^ experimental findings at high concentra- 
tion indicated the relevance of the dipolar interactions^. 
More recently, sensitive echo measurements gave precise 
information on pairs of interacting impurities that oc- 
cur at low doping^. At impurity concentrations of a few 
hundred ppm, the dipolar couplings destroy the coher- 
ent motion of the impurities*. Though a few questions 
would seem settled,^ the role of interaction, especially 
the dynamical aspects, is still controversial. "'^""-'^^ 

At still higher concentration, one obtains mixed crys- 
tals that show glassy behavior with respect to the rota- 
tional motion of the impurities, such as KBr(i_c)(CN)c^'' 
or (MF2)(i_c)(LF3)c, with M=Ca,Sr,Ba. For this latter 



example of mixed fluorite crystals, "'^^ the coupled density 
of states of the impurity degrees of freedom was deduced 
from Raman spectra, at concentrations c = 0.05 — 0.45. 

In the present paper, we study the consequences of in- 
teractions between substitutional defects for the low T 
properties of defect crystals. We are interested in de- 
scribing the effects of interactions on quantized defect 
motion, as well as in analyzing the crossover to glassy 
behavior at sufficiently high defect concentration. Inter- 
actions may be of an electric dipolar nature, or mediated 
by elastic strain fields. 

We consider TV impurities that are randomly dis- 
tributed on a lattice with Nq sites. The defect concentra- 
tion c = N/Nq is usually much smaller than unity; the 
most interesting physics typically occurs in a range from 
c — 10^^ to a few per cent. 

The impurity Hamiltonian comprises a one-particle 
crystal field potential that is identical for each defect, 
and an interaction term that reflects the random config- 
uration on the host lattice. The one-particle potential for 
a given impurity is written most conveniently in terms of 
local coordinates. For Li impurities in a KCl host, these 
coordinates could be chosen as Vi = Ri — where Ri 
is the impurity position at the lattice site ii^. Tunneling 
arises from degenerate minima at off-center positions of 
the crystal field potential G{vi). The local displacement 
\vi\ is significantly smaller than the lattice constant. For 
other systems, the Vi might describe orientational de- 
grees of freedom of a defect molecule. 

In a dipole approximation, the interaction of impuri- 
ties i and j is linear in the local variables Vi and Vj and 
involves a coupling parameter Jij . It may arise from elec- 
tric dipole moments Pj = gv^, or from elastic quadrupole 
moments; in both cases the interaction parameter is 
proportional to the inverse cube of the impurity distance 
ry ^\R°~R% that is - r'^ . 

In our investigation of collective effects we shall, in 
what follows, resort to a simplified description in terms 
of an effective scalar mean-field analysis that consists 
of the following approximations, (i) The dependence of 
Jij on the relative angles of Vi, Vj, and — is 
simplified to two possible signs, (ii) The local coordinates 
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Vi arc replaced by the one-dimensional scalar variables Vi , 
the onc-particlc crystal field potential by G{vi), and the 



dipolar interaction by Ji- 



In this scalar version, the 



impurity Hamiltonian takes the form 



with 



(1) 



(2) 



The distribution Pq of dipolar couplings can be eval- 
uated for spacings rij much larger than the lattice con- 
stant, i.e., for small impurity concentration c -C 1. In a 



continuum approximation, the Jy 
lates into 



Po{Jij) 



1 



2(iVo-l) J' 



scaling trans- 



(3) 



for Jmax/iVo < \Jij\ < Jn 



The factor ^ accomits for 
is 



the two possible signs of Jtj. The upper bound Jn 
determined by the interaction of closest neighbors, 



3g2 



3 • 

min 



(4) 



where we have used the form of the dipole moment 
Pi = qvi. (Note that rmin is of the order of the lattice 
spacing.) Since we have replaced the dipolar interaction 
in three dimensions by the simpler expression JjjVjVj, 
the numerical constant in (4) is to some extent arbitrary. 
Using a different scheme^^ we would obtain 2/3 instead 



of the factor 3/47r. Because of the Ji, 



scaling, 



the lower bound Jmax/-^o gives the coupling of impuri- 
ties whose distance is of the order of the sample size. 

(iii) Lastly, we neglect the distance dependence of the 
couplings and approximate the random-site character of 
the defect problem by a mean field model with randomly 
distributed all-to-all connections. In a further simplifying 
step, the Jy -distribution is taken to be of the Gaussian 
form 




2J2 



with zero mean J^- = and finite second moment. 



(5) 



(6) 



which is determined in such a way that it coincides with 
that of the actual distribution function PQ^Jij), 



j dJijjfjPoiJij) = jyj^-^r 



(7) 



When putting J'^/N = J^^^/Nq and using c = N/Nq, we 
obtain the scaling of the parameter J with defect concen- 
tration c, 



J = \/c Jn 



(8) 



The Gaussian distribution P is roughly constant for | Jij \ 
smaller than J/\fN and vanishes rapidly for higher val- 
ues, whereas the more precise function Pq shows a power 
law behavior. We repeat that, with (8), their second 
moments are identical. 

Within our scalar approximation, we take crystal sym- 
metries into account by demanding that the one-particle 
potential satisfies the condition G{vi) = G{—Vi). More 
specifically we choose G{vi) to exhibit degenerate minima 
at finite displacement Vi = zLa. As the simplest choice of 
an on-site potential with these properties we use 



G{v)=g{v^-ay 



(9) 



The coupling constant g is chosen in such a way that the 
tunnel splitting Aq for impurities moving in an isolated 
double well of the form (9) reproduces experimentally 
observed results for the lowest excitation energy of an 
isolated defect. E.g., in the case of ^Li in a KCl host, one 
would require Aq ~ 1.1 K, given a value of a ~ 0.7 A. 

Due to the simplifications introduced above, the in- 
teraction energy (2) of the defect Hamiltonian has been 
made to resemble that of the SK spin-glass model, "'^^ al- 
beit one with continuous degrees of freedom rather than 
Ising spins. Models of this type have recently been 
proposed as candidates for describing low temperature 
anomalies of glassy and amorphous systems. Indeed, it 
has been shown^^'^* that a frustrated interaction of the 
form considered here is able to produce a potential energy 
landscape comprising an ensemble of single- and double 
well configurations - the latter with a broad distribution 
of asymmetries and, in the translationally invariant case, 
also barrier heights, even if one starts out with single-site 
potentials G{v) which are a-priori of single well form. 

Our model of interacting impurities can be analyzed 
exactly within mean-field theory, i.e. a self-consistent 
representation of the interaction energy as a sum of ef- 
fective single-site potential energies 



Ui^t{{vi})^Y^U,sivi) 



(10) 



becomes exact in the thermodynamic limit. The ensem- 
ble of effective single site potentials Ues{vi) represents 
the potential energy landscape of the system of interact- 
ing defects. It will be found to contain randomly vary- 
ing parameters whose distribution can be computed}'' 
By quantizing the defect motion within the collectively 
determined ensemble of effective single site potentials 
Ucs{vi), one finally obtains a semiclassical description of 
interaction effects on the behavior of quantum impurities. 

We have organized the remainder of our material as 
follows. In Sec. II we analyze the potential energy land- 
scape of the interacting system within mean-field theory. 
The analysis follows a proposal previously advocated for 
glasses. Sec. Ill describes the analytic solution of the 
mean-field theory in the weak coupling regime. Thermo- 
dynamic consequences of interaction effects are explored 
in Sec. IV, and simplifying features of two-state and 
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WKB approximations are discussed in Sec. V. Sec. VI 
is devoted to dynamic effects, in particular to a com- 
putation of the distribution of relaxation rates, and to 
an analysis of the dynamic susceptibility. We discuss 
our results in some detail in Sec. VII, comparing them 
with those of complementary approaches, and with ex- 
periments, and close with a brief summary in Sec. VIII. 



II. MAPPING OUT THE POTENTIAL ENERGY 
SURFACE 

To map out the potential energy surface of the inter- 
acting system, one computes the configurational free 
energy 

MP) = -imr'ln I n'^^iexp[-/3C/pot({^;i})] , (11) 

rising replica theory to average over the ensemble of ran- 
dom Jij matrices, so as to get typical results. The T = 
limit is eventually taken to select one of the (possibly 
many) classical ground-state configurations of the inter- 
acting system. 

As announced above, a mean-field decoupling produces 
a collection of independent single- site potentials Ues{vi) 
with random parameters (which are comparable to ran- 
domly varying local fields in the context of spin-glasses). 
Technically, this decoupling within replica theory can be 
seen as a method for the self-consistent determination of 
the distribution of these random parameters. 

We shall not repeat here the details of such a calcula- 
tion, as they follow standard lines of reasoning^^. One 
obtains /(/3) = lim„_>o /«(/?) for the quenched free en- 
ergy, with 

a,b 

-r ' In / n dv"" exp [ - mes{{v''})] . (12) 

Here 

C^eff(K}) = -^/3^'E««''^"^'' + E^(^") (13) 
a,b a 

is an effective replicated single-site potential, and the or- 
der parameters qat = determined as 
solutions of the fixed point equations 

qab = {v''v^) , 0,6= l,...,n , (14) 

where angular brackets denote a Gibbs average corre- 
sponding to the effective replica potential (13), and where 
it is understood that the limit n — > is eventually to be 
taken. 

We arc, in what follows, going to solve the self- 
consistency equations only within the so-called replica 
symmetric (RS) ansatz for order parameters 



The two order parameters of the RS ansatz must satisfy 
q={{v^))., q = {{vf)., (16) 

in which {. . .)z denotes an average over a zero- mean 

unit-variance Gaussian z while (. . .) without subscript 
is a Gibbs average generated by the effective replica- 
symmetric single-site potential 

Uks{v) = -J^zv - + G{v) . (17) 

with C = f3{q — q). The RS single site potential contains 
a Gaussian random variable z. A sum of single site po- 
tentials randomly drawn from the RS Gaussian ensemble 
(17) constitutes the mean- field representation (10) of the 
potential energy landscape in terms of an ensemble of 
eff'ective single site potentials. Before analyzing this en- 
semble in greater detail, let us note its two most salient 
features. 

First, due to a collective effect mediated by the interac- 
tion, there is a systematic deepening of the double-well 
crystal field potential experienced by the defects. Sec- 
ond, there is an induced distribution of asymmetries ow- 
ing to the linear contribution to [/rs(v). Note that the 
deepening of the double well potential will in particu- 
lar give rise to a renormalization of the tunnel splitting: 
Aq ^ Aq < Aq, reducing the smallest excitation energy 
that occurs in the system of tunneling impurities. How- 
ever, due to the spectrum of asymmetries there will also 
be a spread of excitation energies towards larger values, 
as we shall describe in greater detail below. 

Either of the fixed point equations for qoi q above may 
be replaced by one for C = f3{q — q) 

which turns out to aquire a finite value in the (3 oo- 
limit that is of interest to us here. The RS free energy 
is 

-r' (^^^J dvexp [ - PUrs{v)] ^ ■ (19) 

As the f3 oo limit is taken, Gibbs averages gener- 
ated by Urs{v) are dominated by the value(s) of v which 
minimize U^_s{v), which we denote hy v = v{z) (display- 
ing its dependence on the value of the Gaussian z). An 
immediate consequence is that q = q in this limit, pro- 
vided that q ^ 0. The T = fixed point equations are 
then 

with i) = v{z) minimizing Ujis{v), i.e., to be determined 
as the appropriate solution(s) of 



qaa = q , qab = q , a^b. (15) 



G'{v) = Jy/qz + J^Cv , 



(21) 
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the prime denoting differentiation with respect to v. The 
solution v(z) of (21) is a smooth fimction of z, except at 
z = where v{z) lias a jump-discontinuity, owing to the 
fact that among the solutions of (21) we have to choose 
the one corresponding to the absolute minimum of Z/rs- 
For 2; 7^ then, we obtain 



1.25 



= M 



(22) 



Taking the jump discontinuity of v{z) , hence the S— 
function singularity of its 2;-derivative at 2; = into ac- 
count, we obtain the following form of the fixed point 
equation for C 



C = 



1 



G"{v) - J^C 



Av(0) 



(23) 



where A'()(0) denotes the size of the jump of v at z = 0. 
The T = limit of the free energy, i.e., the internal 
energy u is 



u=\j^C{q + q) + {Uns{v{m. 



(24) 



Owing to symmetry, the fixed point equations always ad- 
mit of a g = solution. In a situation, where G{v) is of 
double-well form (with minima at ±a), however, only 
the q =/= solution is thermodynamically acceptable as 
T ^ 0. Using the fixed point equations, we find that the 
internal energy is given by 



-J^Cq+{G{v{z)))^ 



(25) 



in this limit. 



III. THE LOW-CONCENTRATION OR 
WEAK-COUPLING REGIME 

The weak-coupling (low-concentration) limit is for- 
mally defined by the inequality Ja'^/ga^ <C 1. It states 
that typical interaction energies are much smaller than 
the bare classical barrier height of the on-site poten- 
tial G{v). In this limit, one has approximately v{z) ~ 
asgn(2;) for moderate values of the Gaussian z, hence 
q = q c^i a?' and C ~ J~^^J2/'k. Here ±5 denotes the 
coordinates of the minima of U-rs{v) at 2 = 0, that is 



a = a 



1/2 



(26) 



In the weak coupling limit, therefore, the effective single 
site potential experienced by the defects is of the form 



-Jaz V 



J 



2n 



G{v) , 



(27) 



exhibiting both, the systematic deepening of the double 
well potential experienced by the defects, and the induced 
distribution of asymmetries due to the linear contribution 
to J/rs {v) mentioned above. Both scale linearly with J, 



10° 



10° 



10"* 



10"' 



FIG. 1. Test of low-concentration approximation. Upper 
set of curves: q, lower set of curves: JC, given as functions of 
the defect concentration c. Deished lines give the results of the 
low-concentration approximation, full curves, the numerical 
solutions of the full fixed point equations. 

hence with the square root of the defect concentration in 
the low concentration regime (see (8)). 

As Figure 1 shows, the results of the low concentra- 
tion approximation agree fairly well with those of a full 
numerical solution of the fixed point equations for all 
concentrations of interest. 



IV. THERMODYNAMICS AT LOW 
TEMPERATURES 

The contribution of the defects to the thermodynamics 

of the system at low temperatures is dominated by quan- 
tum effects. That is, one has to consider the quantized 
motion of the defects in the effective single site potentials 
given by (17) or by their weak-coupling approximations 
(27). The first task therefore consists in determining the 
energy levels and the corresponding eigenstates ipn 
for the ensemble of effective single site Hamiltonians 



2m 



(28) 



This done, one may proceed to compute densities of state. 
Normalized with respect to the total number A^o of lattice 
sites, this gives 



Pn{E) = C {d{E - En) 



(29) 



where En = En — Eq. The defect contribution to the 
specific heat in the same normalization gives 



c = ckBPmnls)-{nesf)^ 

in which now 

T^(...)exp(-/3Weff) 



(...) 



Trexp(-/?Weff) 



(30) 



(31) 



Another quantity of interest is the static (dielectric) sus- 
ceptibility, originating from the interaction of the defect- 
dipoles with an external field £, Hext = —q£ Vi. The 
static Kubo-formula gives 
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The v„m = {'>Pn,'vipm) dcnotc matrix elements of the 
position operator between the various eigenstates of the 
effective Hamiltonian Hes, and Z^e is the canonical par- 
tition sum generated by it. 



10"' 

10-= 




1 10 100 1000 

E[K1 



FIG. 2. Density of states for the five lowest lying excitation 
energies at c = 10 ppm 

As shown in Figure 2, the lowest band of excitation en- 
ergies, originating from tunneling-excitations in the en- 
semble of asymmetric double-well potentials extends to 
much lower energies than the other excitations. The lat- 
ter are related to harmonic excitations about the two 
minima of G with energies Uwo only beyond 100 K. The 
pair of bands with peaks near 99 K and 125 K again corre- 
sponds to a pair of states mixed due to tunneling between 
the wells. The other two bands (with peaks near 192 K 
and 260 K correspond to states with energies above the 
barrier. 

Tunneling excitations within the space spanned by the 
oscillator ground states in the two wells will dominate 
the low temperature physics. This is nicely seen in the 




E[KI 

FIG. 3. Density of states for the lowest band of tunneling 
excitations at c = 1, 10, 100, and 1000 ppm (bottom to top). 



specific heat data exhibited in Figure 4. Notice that the 
Schottky peak at low concentration is modified through 
interaction effects. As c is increased beyond 1000 ppm, 
a range of temperatures develops where the specific heat 
starts to show a linear temperature dependence much 
like in glasses. This is due to the fact that the density 
of states corresponding to the (lowest band of) tunneling 
excitations shown in Figure 2 becomes nearly constant 
in the energy range Aq < E < ^/cJ^i^^a^, which covers 
a large range of energies, as the concentration is suffi- 
ciently increased; see Figure 3 which also exhibits the 
renormalization of the tunneling matrix element towards 
lower energies due to interaction effects. 

10-' = — — — — 



10-^ 

c 

10-' 




T[K] 

FIG. 4. Specific heat as a function of temperature for c = 
1, 10, 100, and 1000 ppm (bottom to top). The influence of 
higher order excitations is seen only at temperatures above 
10 K. The dashed curve is a Schottky-peak. 

In the static susceptibility plotted in Figure 5, no sig- 
nificant contribution of higher order excitations is de- 
tectable (on the scale of the figure) at all. even up to tem- 
peratures as high as 60 K. The same quantity evaluated 
in a two-state approximation would be indistinguishable 
from what is shown here. Note that the concentration 
dependence at low temperatures is proportional to c at 
low concentrations, but crosses over to a -^/c behavior at 
larger concentrations. We shall return to this in greater 
detail later on. 



V. TWO-STATE APPROXIMATION 

At low temperatures, kBT <C huiQ, only the two lowest 
lying states are significantly populated. For a sufficiently 
high barrier they can be constructed in terms of a basis 
of pocket states \L) and \R) that are localized about the 
potential minima at ±a of the two wells. 

The corresponding Hamilton matrix reads 

where the off-diagonal elements arc given by the renor- 
malized tunnel energy. The diagonal entries account 
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for an asymmetry energy between the two wells, A = 
pRsi-a) - Unsia)]; with (27) we find 



2Ja'^z = Az . 



(34) 



The two-state Hamiltonian is easily diagonalized. Its 
eigenvalues are given hy ±^E with 



+ A2; 



(35) 



the corresponding eigenstates |±) may be written in 
terms of a mixing angle tant^ = Aq/A, 



1+) = cos (j)\L) +sin^|i?), 
|— ) = sm(j)\L) — cos(j)\R). 



(36) 



Note that the asymmetry energy A is linear in the 
Gaussian variable z; its distribution law is thus given by 



Pa(A) 



1 



27rA 



exp 



2A2 



(37) 



As a consequence, both the two-state energy splitting E 
and the mixing angle (j) are spread over a range deter- 
mined by this distribution function. 

The resulting density of states for the lowest band of 
excitation energies {p{E) = pi (E) in the notation of the 
previous Section) 



p{E) 



6 E- 



+ A2 



(38) 



is easily calculated. As is obvious from (35), we have 
p{E) =0 fov E < Aq] at larger energies it is determined 
by the distribution of asymmetries, so 



p{E) 



2E 



7^ 



Pa \ E^ 



(39) 



for E > Ao, and shows a square root singularity for E — !■ 
Ao; the pre- factor 2 is due the restriction E > 0. 

At very low concentration, the width of Pa (A) is much 
smaller than the tunnel energy Aq, and p{E) is sharply 



0.00004 



0.00003 



X 0.00002 



T[K] 



FIG. 5. Static susceptibility as a function of temperature 
for the same concentrations as in Figure 3. 



peaked about Aq. As soon as the typical interaction 
energy A:^2J a? exceeds the tunnel energy Aq, p{E) is 
well approximated by a Gaussian above Ao. It remains 
to compute the renormalized tunnel energy Aq. 

With respect to the computation of Aq, we have looked 
at two variants. In the first, both Ag and its renor- 
malized value Aq are computed exactly by solving the 
Schrodinger equation of the impurity moving in the ap- 
propriate crystal field potential, i.e., G{v) for the former, 
and the renormalized potential G{v) ~ —-j=v'^ + G(v) 
for the latter. As mentioned above, the first computa- 
tion is actually used to fix the coupling g in (9) so as 
to reproduce the experimentally observed value - e.g., 
Aq I.IK for ''Li in KCl. With parameters appropriate 
for the KCh^Li example, this requires a potential with 
bare barrier height go?' ~ 174K. 

In the second variant, both Aq and its renormalized 
value Ao arc computed within a WKB approximation. 
The bare tunneling matrix element is in this setting usu- 
ally parameterized as 



with 



Aq = TujJq exp(— A) 



X = d^2mVB/h^ 



(40) 



(41) 



and parameters derived from the bare on-site potential 

G and the mass m of th e tunneling particle — the fre- 
quency Wo = v^Sga^/m of harmonic oscillations about 
the minima of G, the distance d = 2a between the min- 
ima, and the height Vb = ga"^ — hwo/2 of the classical 
barrier above ground state in the two wells. Once more, 
this first calculation would be used to fix the coupling 
constant g of the bare potential (9). Within this approx- 
imation, the classical barrier required to reproduce the 
bare timncling energy reported for KCh^Li is ga^ ^ 84 K, 
i.e. only approximately half the value obtained from the 
numerically exact analysis. The renormalized tunneling 
matrix clement is computed in the same way, except that 
parameters are computed from the renormalized on-site 
potential. This gives 



(Do = oJoa/a 



Vb = Vb + - 
d = 2a 




(42) 



to lowest order in small parameters (with a /a ~ 1 -|- 
J/iA^ga"). 

The renormalized tunnel energy Aq = luboe~^, with 



\ = d^2mVB/n^ 
may thus be expressed as 



Aq = AoVl + ee" 



(43) 
(44) 



Here 
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FIG. 6. The ratio Ao/Ao as a function of concentration. 
Exact evaluation (full line), WKB approximation (dashed 
line). 



2 Jo? 



and 



^1 



huJn 



(45) 



(46) 



the latter quantity being evaluated in the weak coupling 
approximation, Ja^ ^ ga^ which is always appropriate 
for the case at hand. If, moreover, the barrier is assumed 
to be high, ?iuo <^Vb, so that there are many oscillator 
levels between the potential minima at w ~ ±a and the 
top of the barrier at t; = 0, the expression simplifies to 



= a| J-:^ :^ 3./-%^ , (47) 
2 8\ TV ga^ y n JkJo ^ ' 



in which corrections involving powers of Ja? /ga^ and 
huio/ga'^ have been neglected. There is no restriction 
on the ratio Ja^/huio- It turns out, however, that in the 
KCl:Li case the ratio hujo/ga'^ is roughly 1.2, thus not 
small, so that the corresponding simplifications are not 
available. Figure 6 compares the ratio Aq/Aq evaluated 
numerically and via the WKB approximation (40), (41) 
and (44)- (46). The renormalization of the timnel energy 
becomes noticeable when the concentration exceeds 100 
ppm, a concentration at which the typical asymmetry 
A becomes comparable with the vibrational energy fkoo 
in the two wells; the renormalization effect is, however, 
overestimated within the WKB approximation. 

In the two-state approximation, specific heat and static 
susceptibility are given by 



C = ckB j dEp{E) 



cosh^{f3E/2)' 



(48) 



and 



X = cl3{qa) 



\cosh^{l3E/2) 



+2/cbT-^ tanh(/3^/2)\ , (49) 

/ z 

respectively. The results are valid as long as tempera- 
ture is much smaller than the librational energy huo. In 
(49) we have also introduced the usual approximate rep- 
resentation V = CLGz of the position operator in a basis 
of pocket states. This latter approximation is the largest 
source of errors in the expression (49) for the susceptibil- 
ity in that it overestimates matrix elements by roughly 
10%. It turns out tat this has a noticeable effect only on 
an overall pre-factor (at least for isolated or weakly cou- 
pled defects), but not on the temperature dependence. 



VI. DYNAMICS 

A. Phonon damping 

The interaction of an impurity at site i with elastic 
waves is described by the coupling potential 



(50) 



with the elastic deformation potential 7, the reduced co- 
ordinate Vi/a, and the phonon strain field 



2Vquj. 



qs 



-qi {bqs -b^s) , 



(51) 



where g is the mass density of the host crystal, and q 
labels the wave vector of three acoustic phonon branches 
s. There is ample evidence that the defect-phonon cou- 
pling is weak; therefore it may be treated in first Born 
approximation. 

All dynamic information may be obtained from the 
two-time correlation function 



G{t) = {l/2{v')){v{0)v{t) +vit)v{0)), 



(52) 



with the normalization condition G(0) = 1. At low tem- 
peratures, we may use the two-state approximation. We 
shall here restrict ourselves to the simplified represen- 
tation V = aaz of the position operator in the basis of 
pocket states as introduced above; it is sufficiently pre- 
cise to give qualitatively reliable results. 

The theory of a weakly damped two-state system has 
been derived in many places; here we merely quote the 
result for the correlation function, 



G{t) = -^cos{Et/n)e-^ 



Ft 



+|J((l-Q)e-" + Q) 



(53) 



with Q = tanh(/3iS/2)^ and the one-phonon damping 
rate 



T=^^AlEcoth{f3E/2). 



(54) 
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Here Vs denotes the sound velocity. G{t) determines the 
linear response of a two-state system with dipole moment 
p = qaaz, tunnel energy Aq and asymmetry A. There are 
several interesting issues arising from the above model. 

First, the Gaussian distribution for the bias A leads 
to a distribution of resonance energies E and relaxation 
rates F. 

Second, because of the variation of the typical asym- 
metry A with concentration c, the relaxation behavior 
of the impurities changes significantly with increasing c. 
Prom (53) it is immediately clear that there is no sig- 
nificant relaxation contribution for A <C Aq, i.e., at low 
concentration. In the opposite case A Aq, the os- 
cillatory (resonant) part of G{t) is insignificant, and the 
relaxation term governs the impurity dynamics. 

Third, the temperature factor of the relaxation con- 
tribution, 1 — Q = cosh.{PE/2)~'^, vanishes at very low 
temperatures ksT <^ E, whereas it tends towards unity 
at higher T and large concentration. 



B. Rate distribution 

Recent measurements of the dielectric constant of 
KCL:Li revealed relaxational motion of the lithium im- 
purities over several decades in the kHz range. ^ Such a 
broad relaxation spectrum is characteristic for disordered 
system in general. Therefore we discuss in some detail 
the distribution of relaxation rates 



(55) 



In our two-state model there is only one variable, A = 
zK, which is entirely determined by the Gaussian distri- 
bution of z and the constant A; cf. (35). 

We give explicitly the rate distribution at zero tem- 
perature, because of its simplicity and since it shows the 
essential features. At T = 0, Eq. (54) gives 



with 







T=l^E 
Ao 



and Fn = 



1 37^ 



rAL 



(56) 



(57) 



Hence, except for scales, the rate distribution is at the 
presently chosen level of description basically equivalent 
to the density of states in the two-level approximation. 
Substituting the Gaussian variable z appearing in the 
energy E, we find 



1 Ao 



2F 



27r AF 



Y2 



F2 
^ 



exp 



^0 ^p2 



2A2F2 



(58) 




FIG. 7. Rate distribution at T = for c = 1, 10, 100 and 
1000 ppm. 



Here, Fq is the minimum rate at finite doping. Ac- 
cording to (57), it is reduced by a factor Aq/Aq^ as com- 
pared to the rate in the dilute limit, Fq. Like the density 
of states, -Pr(F) contains a Gaussian factor and a square 
root singularity at the minimum value Fq. The latter is 
of little relevance, and the rate distribution is governed 
the exponential in (58). We point out the most salient 
features of the rate distribution. 

(i) For low impurity concentration, we have Aq ^ A, 
and ^'r(F) is a sharply peaked function with the lower 
cut-off Fq. In the dilute limit c ^ 0, the distribution 
tends towards a delta function at Fq. However, already 
at concentrations as low as 1 ppm, the rate distribution 
has a visible tail towards higher rates due to the presence 
of asymmetries (see Figure 7). 

(ii) In the opposite case of high doping, the maximum 
asymmetry exceeds by far the reduced tunnel energy, 
Ao <i; A. As a consequence, the Gaussian factor in (58) 
leads to a wide distribution that is almost constant be- 
tween the lower bound Fq and the effective upper cut-off 
(A/Ao)fo. 

(iii) The lower bound of the distribution, Fq, depends 
on the impurity concentration. In the dilute case, the 
tunnel energy is given by the bare value Aq. Yet at finite 
concentration, the renormalization of the tunnel energy 
reduces the rate by the factor Aq/Aq^. At c = 0.01, 
the minimum rate is by roughly one order of magnitude 
smaller than in the dilute case. (The WKB approxima- 
tion predicts more than two orders of magnitude; see 
Figure 6). Thus interaction leads to both slow and very 
fast relaxation, as compared to the low-doping case. 

Figure 7 shows the evolution from low to high doping. 
For concentrations up to 100 ppm, the curves have the 
same shape as those for the density of states (as predicted 
within the two- level approximation using the pocket state 
approximation for the position operator. The curve for 
1000 ppm develops a different shape at high rates, be- 
cause the pocket state representation of the position op- 
erator becomes less and less precise at large asymmetries 
(compare Figures 7 and 3). 



for F > fo. 
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C. Dynamic susceptibility 

Experimental investigations of the two-state dynam- 
ics in terms of clastic or dielectric response function in- 
volve the dynamical susceptibility x(<^) = x'('^)+* x"!'^); 
whose spectral function 



X"(^) 



: c - tanh(/3;ia)/2) (G"(w)), 



(59) 



is related to the average of the motional spectrum G"{u}), 
i.e. of the Fourier transform of (53). (The additional 
prefactor (qa.)'^ appearing in the susceptibility (49) and 
(32) is due to the fact that in Sees. Ill and IV we have 
considered an external field coupling to the dipole oper- 
ator qv_ra,thev than to a normalized position operator 
v/^y (f 2) whose correlator is considered in (52)). 

Here we give explicitly the real part that describes the 
sound velocity or the reactive part of the dielectric func- 
tion. For relevant frequencies fku <^ Aq (and relaxation 
rates at typical concentrations satisfying F^ <C Aq) one 
has 



X'(c^) 



X^tanh(/3i?,/2) 



A? 



(60) 



within the two-state approximation. There are two con- 
tributions of different origin to the susceptibility. The 
first or "resonant" part is dominant in the dilute limit, 
where the typical asymmetry is significantly smaller than 
the tunnel energy, A <C Aq. On the other hand, the 
second contribution prevails at strong doping, where 

A > Aq. 

In Fig. 8 we plot x'(a;) as a function of temperature 
for various impurity concentrations c. At low doping, 
the susceptibility is of the van Vleck type and hardly de- 
pends on the external frequency. The temperature vari- 
ation is given by the occupation difference of the two 
levels, tanh(Ao/A:T); accordingly x'{^) is constant for 
kT < Aq. At larger doping, pronounced relaxational 
contributions occur, with peak-positions strongly depen- 
dent on frequency. 

The zero-temperature value of x' is a convenient mea- 
sure of the relevance of interaction effects. The temper- 
ature factor in the relaxation contribution vanishes at 
T = 0, and one has 



x' = c^((l + zW/Ag)-3/2), 



(61) 



The average over z can be written in terms of the con- 
fluent hyper-geometric function U{a, b, x), 



X =c- 



-.U 



1 1 

2'"' 2^;2 



(62) 



where we have defined the dimensionless coupling param- 
eter f = (A/Ao), 



10 




0.1 



1 

T[K] 



10 



FIG. 8. Real part of the susceptibility as a function of tem- 
perature at tj = 10 kHz for c= 10, 100, 1000 and 10000 ppm 
(bottom to top). For the c = 1000 ppm curve, the resonant 
(dashed) and relaxational (dot-dashed) contribution are sep- 
arately exhibited as well. 



Vc2 



~ Vc2- 



(63) 



For weak doping, the asymmetry A is much smaller than 
the tunnel energy; with < 1 and U{l/2,0,l/2iy'^) « 
•\/2j^ we recover the obvious result x = c/Aq. In the 
opposite case of strong doping we have i/ ^ 1 and 
f/(l/2,0,l/!/2) = = 1.128... From (62) we thus ob- 
tain in the limiting cases 



X 



c/Ao for c < Co 

c/A = U(x)\/cco/Ao forc»co 



(64) 



The cross-over occurs at a concentration cq where the 
asymmetry attains the value of the tunnel energy, A w 
Aq, or z/ Ri 1. Eq. (62) yields 



Co = 



2 T 



(65) 



Because of the concentration dependence of Aq this is an 
implicit equation for cq. With parameters appropriate for 
''KCLLi as used before, one finds cq = O{10~^); compare 
Figure 9. 

The relaxation contribution to (60) strongly depends 

on the frequency lu and the rates Fj. Relaxation is most 
efficient where the rate of thermal two-level systems are 
close to the external frequency u; lowering the frequency 
shifts the peak to lower temperatures. The relaxation 
peak is significantly broadened by the rate distribution 
(58). 



VII. DISCUSSION 

A. Cross-over to relcLxation 

Broad spectra of energies and relaxation rates are char- 
acteristic for any disordered system. In the case of quan- 
tum impurities with dipolar interactions the broadening 
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is tuned by the parameter which measures the strength 
of the interaction in units of the tunnel energy and in- 
creases with the square root of the impurity concentra- 
tion c; cf. (63). 

At low doping the motional spectrum is peaked about 
the bare tunnel frequency Ao/7i. The relaxation rates 
that are close to Fq as defined in (57) are of little im- 
portance, since the weight of the zero-frequency feature 
of the susceptibility is small. This changes as the impu- 
rity concentration reaches the value Cq, i.e. as v tends 
towards unity. Then the average asymmetry energy is 
comparable to the timnel energy; as a consequence, a 
strong relaxation peak emerges, and both the density of 
states p{E) and the rate distribution ft(r) broaden; cf. 
Figs. 3 and 7. 

The static part of the susceptibility at T = provides 
the clearest signature of this cross-over to relaxation. In 
Fig. 9 we compare the result of the present work (62) 
with that obtained previously by one of us in a Mori 
projection scheme^^, 



X =Xo- 



( yrT7?-M)^ 



(66) 



f 



0.01 



0.001 



0.0001 




where the dimensionless coupling constant 



0.0001 



FIG. 9. Zero temperature susceptibility as a function of the 
concentration. The solid line gives the result of the present 
work, the dashed line that of Rcf. 13, i.e. Eq. (66), and the 
dotted line indicates non-interacting impurities. Full symbols 
indicate experimental data on KCli^Li and KCb'^Li [ 20, 9, 6], 
open symbols on KC1:0H [ 6, 21, 22], LiF:OH [ 21], RbCl:OH 
[ 6], NaChOH [ 21]. (In order to obtain comparable quantities 
for these various systems, we plot the rescaled susceptibility 
A^x/Xo as a function of the dimensionless quantity ji that is 
linear in c.) 



H = c2- 



(67) 



may be considered as the rescaled concentration. (The 

extra factor 2 in (67) as compared to the definition in 
Rcf. 13 is due to the fact that our convention for the cou- 
pling between impurities differs from that used in Rcf. 13 
by a factor 1/2). The susceptibility in the absence of in- 
teraction, xo = c/Ao, varies linearly with the concentra- 
tion; in Fig. 9 it is indicated as a dotted line. Moreover 
we have indicated experimental data for various doped 
alkali halides. 

At very low doping, both expressions (62) and (66) 
are linear in the concentration, as expected for non- 
interacting impurities, whereas at higher impurity den- 
sity the dipolar interactions of the tunnel systems 
strongly reduce the susceptibility. We discuss the differ- 
ences of the theoretical results obtained from the present 
mean-field model (62) and the Mori projection scheme 
(66). 

(i) Both parameters and v are proportional to the 
ratio of the maximum interaction JmaxCi^ and the tunnel 
energy. In the present approach we found a square root 
dependence v oc whereas projection scheme yields a 
linear law jj, (x c. (Note that in the expression for u we 
have neglected the weak c-dependence of Aq.) 

(ii) As a consequence or these different powers, the 
cross-over to relaxation occurs at much higher concen- 
tration in (66), « 1 or Co = Ao/JmaxO^; the present 
mean- field theory yields, at f w 1, the square of this 
quantity; cf. (65). 

(iii) In the strong-doping limit, the susceptibility in 
(62) varies as x' \fc-\ it still increases with concentra- 
tion, whereas (66) results in a decrease, x' c~^. 



These discrepancies are easily traced back to the ba- 
sic assumptions of the models. The present mean-field 
approach relies on random couplings to all impurities. 
The characteristic interaction J is given by the second 
moment of the Gaussian distribution and thus is propor- 
tional to the square root of the concentration; cf. (8). 

In the previous Mori approach^"^, the characteristic in- 
teraction is determined by a few nearest neighbors of a 
given impurity. (Because of the r~^-dependence on dis- 
tance of the dipolar interaction and because of its ran- 
dom sign, neighbors at larger distances contribute little 
to field at a given site.) The r~^-law results in the linear 
variation of J with concentration. 

It should be noted that the square root law (8) can be 
modified by a change of the mean-field model in a man- 
ner that takes the dominant effect of the interactions of 
a defect with a few nearest neighbors into account. In 
a mean-field context this might be achieved by replac- 
ing the fully connected random interaction matrix by a 
sparse random matrix such as in the Viana-Bray modeP^, 
in which typical interaction scales may be chosen to scale 
linearly with the concentration. In the details, however, 
such a diluted mean-field model is much more compli- 
cated, requiring the introduction of infinitely many order 
parameters. 



B. Comparison with experiment 

The most significant result of the present work consists 
in the cross-over at c = cq from coherent one-particle 

motion to relaxation driven by the dipolar interaction. 
(The pre-factor (A/i?)^ of the relaxation feature in the 
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susceptibility (60) vanishes for zero interaction, whereas 
it tends towards unity in the strong-couphng hmit.) 

Such a cross-over and the emergence of a relaxation 
peak in the motional spectrum have been observed for 
various impurity systems such as KCl:Li and NaChOH.^ 
As discussed above and shown in Fig. 9, the concentra- 
tion dependence of the susceptibility in the strong-doping 
limit obtained in the present work, x' cx \/c, differs from 
the previous result x' c^^ (Ref. 13). The data for the 
highest concentrations (/U > 0.1) clearly show the rele- 
vance of the dipolar interactions; the measured values 
are orders of magnitude smaller than expected for non- 
interacting impurities. 

In spite of their marked differences of the theoretical 
laws (62) and (66), the available data do not permit us 
a definite statement on the power law x oc c". Though 
(66) provides a better fit at intermediate concentrations, 
0.001 < ji < 1, the data do not settle the power law at 
high concentrations. (Unfortunately, presently available 
crystals do not satisfy the strong doping criterion /i ^ 1.) 
Moreover, it is rather a difficult matter to properly sepa- 
rate the relaxation contribution for strongly doped crys- 
tals; thus the data points at the highest concentrations 
may well be too large. Further experiments at high den- 
sities would be most desirable, in view of the above dis- 
crepancy between a = 5 in the present approach and 
a = —2 from Ref. 13. 

Wc now turn to the rate distribution plotted in Fig. 
7, that shows a strong broadening with rising concen- 
tration. Indeed, the relaxation spectra observed at high 
concentrations become very broad and develop a signif- 
icant low- frequency wing^. On the other hand, specific 
heat measurements give clear evidence for a broadened 
density of states, in qualitative agreement with Fig. 3. 
(Cf. Ref. 13 and original literature cited therein.) 

Wc close the discussion of alkali halidcs with a re- 
mark on the temperature dependence of the relaxation 
amplitude in the susceptibility (60). Since our mean- 
field theory reduces to a set of effective two-state sys- 
tems, the relaxation feature shows the well-known fac- 
tor cosh(i/3_E)~^ and thus vanishes exponentially in the 
limit T ^ 0, whereas the rate F tends towards a con- 
stant. This behavior is characteristic for a local degree 
of freedom. The Mori projection approach, on the other 
hand, results in a constant relaxation amplitude and a 
rate that decreases at low T, thus showing the generic 
behavior of collective relaxation. Again, the experimen- 
tal situation would seem not entirely conclusive, though 
the data on KCl:Li of Ref. 8 would suggest that both be- 
haviors are present, i.e. the relaxation spectrum would 
comprise one-particle and collective contributions. 

Finally we discuss recent Raman light scattering ex- 
periments on mixed fluorite crystals (MF2)(i_c)(LF3)c, 
with M=Ca,Sr,Ba.^^ The reported data give strong ev- 
idence that the coupled density of states D{E) of the 
two-level systems increases linearly with c in the range 
c = 0.05 - 0.5, i.e. D{E) = cDq{E); it would seem that 
the shape function Do{E) does not change in this range. 

We consider very likely that the broad distribution of 
two-level energies observed in Ref. 15 arises from elastic 



coupling of LF3 tunneling states. In our model the TLS 
coupled density of states of is given by 

D{E) = ^p{E). (68) 

With the expression (39) for p{E) one easily finds that 
D{E) varies over a wide range as E~^ and increases with 
c. Though these dependencies on E and c qualitatively 
agree with the data of Ref. 15, this experiment certainly 
requires a more careful study of the high-doping case. 
Possible discrepancies with the present description for 
this particular case arise from the manner in which defect 
crystals can, or cannot become glassy at high doping, an 
issue to which we now turn. 



C. Crossover to glassy behavior 

The appearance in our theory of a broad spectrum of 
excitation energies and a correspondingly broad spec- 
trum of relaxation rates due to interaction effects at 
high doping is reminiscent of the physics characteristic 
of structural glasses. In this sense, our mean-field ap- 
proach appears to describe a crossover to glassy physics 
at high doping. Indeed, one of the original motivations 
for studying defect crystals has been that they consti- 
tute systems in which - unlike in glasses - the nature of 
the tunneling excitations is well understood, and which 
at the same time offer the possibility of approaching a 
glassy limit by increasing the defect concentration (for 
a recent review, see Ref. 27). It should, however, be 
noted that the type of solely interaction-mediated glassi- 
ness described in the present paper is different from that 
believed to describe glasses proper in one essential aspect: 
while tunneling systems do occur with a broad distribu- 
tion of asymmetries, there is - unlike in glasses - no cor- 
responding randomness in barrier heights. This feature 
entails^* for instance that the typical glassy plateau of 
the internal friction as a function of temperature would 
be absent in systems to which the present theory applies. 
On the other hand, such plateaus are known develop in 
certain defect crystals such as (CaF2)(i_c)(LF3)c, but 
they require defect concentrations exceeding the 10% 
range. There are basically two ways to create the re- 
quired randomness also in the barrier heights. First, 
it can be solely interaction-mediated, if interactions are 
translationally invariant. This mechanism is, however, 
not available for a system of defects embedded in a crys- 
talline host. The other possibility is to allow the crystal 
field potentials G{vi) to vary randomly from defect-site to 
defect-site. For highly doped defect crystals, this would 
seem like a realistic, and indeed expected feature. Un- 
fortunately, though, there are no good theoretical models 
around to predict the associated kind of randomness and, 
in particular, its variation with defect concentration. In 
this sense, glassy defect crystals are apparently not much 
simpler systems than structural glasses proper. 
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D. Comparison with classical reorientation model 

The dissipation rate (54) describes jumps between two 

quantum levels; it may be decomposed as F = + Fj, 
where Ff accounts for thermally activated jumps from 
the ground state upwards and F^ for the reverse process. 
Detailed balance requires F| = e~^^Ti; thus we obtain 
with E ^ kT the thermally activated rate 

F| = foe"'^^, (69) 

where the activation energy is, in the strong-coupling 
limit, given by the dipolar interaction, E = 2Ja?. 
Such a rate, with an appropriate distribution of bar- 
rier heights, has been used in previous work on a clas- 
sical reorientation model for quadrupolar glasses such as 
KBr(i_c)(CN)c^'*'^'*; in the classical picture the reorien- 
tation of the cyanide molecules (i.e. rotation by tt) re- 
quires to overcome their quadrupole-quadrupole interac- 
tion that corresponds to our 2Jo^. Thus the present 
quantum mechanical calculation yields in the classical 
limit the proper temperature dependence of the rates. 
(The present rate distribution is, however, much simpler 
than that derived in Ref. 24.) 

E. Assumptions, Approximations and Potential 
Extensions 

Here we briefly summarize once more the main as- 
sumptions and approximations madc^ in this paper, try- 
ing to assess their quality, and we indicate possibilities 
for further improvements. 

(i) The 3D local coordinates of the impurities have 
been replaced by a scalar variable, and the dipolar cou- 
plings by scalar couplings with a random sign. In the 
analysis of dynamical effects we have introduced a fur- 
ther simplification by truncating the Hilbert space of the 
single-particle effective Hamiltonian to the ground state 
doublet. Whereas the neglect of two additional spatial 
coordinates should be of little relevance, in that the sim- 
plified model still brings out the salient features of the 
problem of interacting tunneling impurities, the trunca- 
tion of the one-particle Hilbert space certainly breaks 
down as soon as the thermal energy or the mean dipo- 
lar interaction exceed the vibrational energy within one 
well of the impurity crystal potential. Yet this is not a 
severe constraint for our model since the most interest- 
ing physics occurs at low temperatures and for dipolar 
couplings comparable to the tunnel energy. 

(ii) Interactions between tunneling impurities are ana- 
lyzed classically at the mean field level. This by itself is 
expected to give a reliable picture for much of the physics, 
as long as we are not interested in coherent quantum mo- 
tion of several partic;les (see below) or critical behaviour 
associated with phase transitions. 

(iii) In actually performing the mean-field analysis we 
have replaced the distribution (3) by the Gaussian (5) 
with the same variance. As discussed in Sect. VII- 
A above, this is certainly a most serious approxima- 
tion, which is dictated solely by the demand for analytic 



tractability, and which can not be justified on physical 
grounds. There would be basically two possibilities to im- 
prove on this within random-bond modeling and replica 
techniques as used in the present paper. One is to con- 
sider sparse connectivity models^^, the other to follow 
a cavity approach proposed by Cizeau and Bouchaud^^. 
However, the former gives rise to an infinite set of order 
parameters and is quantitatively more or less untractable, 
the latter uses the two-state nature of Ising spins in ways 
so essential to the analysis that we have seen no way to 
apply it to continuous degrees of freedom as they occur 
in the present problem. Whether techniques invented for 
models with site-randomness will offer a way out of this 
dilemma, is currently under study. 

Note that the uncertainties with respect to the con- 
centration dependence alluded to above are quite general 
and arise in related approaches to, e.g., quantum spin 
sytems^^. 

A possible and uncontrolled ad-hoc modification of the 

theory which would give a c^/"^ behavior of the T = 
susceptibility at large doping instead of the ^/c behavior, 
and which would therefore fit the high doping behaviour 
of this quantity rather well, is briefly mentioned in our 
summary below. 

(iv) The approximation of replica symmetry used in 
the quantitative analysis of collective behaviour is not ex- 
pected to be a serious source of errors. We have checked 
elsewhere^^'^^ for related systems that replica symmetry 
breaking effects are small e.g. in the low-temperature 
specific heat, but also in distributions of parameters char- 
acterizing double well-potentials. 

(vi) The mean-field model itself constitutes an approx- 
imation that is by no means innocuous for all questions 
one might wish to look at. By mapping the pair inter- 
action on an effective one-particle potential, we discard 
all effects of coherent motion of the impurities. As a 
straightforward consequence, the relaxation amplitude of 
the susceptibility (60) vanishes in the limit of zero tem- 
perature. 

(vii) Finally, we briefiy mention minor approximations 
that are of little significance with respect to the basic 
features of our model. In order to obtain simple laws 
in terms of the concentrations we have truncated various 
powers series in Sect. VI, even if higher-order corrections 
are not always negligible. 

VIII. SUMMARY 

We have studied interacting quantum impurities in 

terms of a mean-field model with a Gaussian distribu- 
tion for the couplings. We briefly summarize our main 
results. 

(i) The dipolar interaction leads to a reduction of the 
tunneling amplitude and to a wide distribution of the 
asymmetry energy. As a consequence, the density of 
states is smeared out to both to smaller and higher values 
as compared to the unperturbed tunnel energy Aq. Sim- 
ilarly, in the strong-coupling limit the relaxation spec- 
trum covers several orders of magnitude. These features 
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account for the collective nature of the underlying re- 
laxation process; they are in qualitative agreement with 
experiments on lithium doped potassium chloride.^ 

(ii) As a clear signature of interaction effects we con- 
sider how the zero- frequency susceptibility at low T varies 
with the impurity concentration c. The present result, 
X oc y/c, strongly deviates from the law x oc found 
previously in a different approach. Though more recent 
experiments on KCl:Li would seem to speak in favor of 
this latter result,^ the data presently available in the 
strong-doping regime /it > 1 do not unambiguously an- 
swer this question. Note, however, that the y/c scaling is 
clearly a consequence of the 'naive' mean-field assump- 
tion of identically distributed all- to-all interactions. It is 
likely to be modified in more realistic approaches. For in- 
stance, by introducing an ad-hoc c scaling of our interac- 
tion parameters in an otherwise structurally unmodified 
mean-field approach, one would obtain a c^/^ behavior 
of the T = susceptibility at large doping instead of the 
■y/c behavior. 

(iii) Since the above square root behavior results from 
a generic feature of the Gaussian mean-field model, an 
experimental test of this law would be interesting for a 
wider class of models. Mean-field models with random 
couplings similar to that considered in this paper are fre- 
quently studied in view of quantum spin glasses and other 
disorder quantum systems.^® 

(iv) Our theory describes a crossover to glassy behavior 
in the (restricted) sense that broad and rather fiat dis- 
tributions of excitation energies and relaxation times de- 
velop at high doping. We have argued that true glassincss 
will develop in defect crystals only, when mechanisms are 
invoked which are beyond those originating from defect 
interactions. As a consequence, true glassy defect crys- 
tals are almost as difficult to grasp theoretically as truly 
amorphous systems. 

ACKNOWLEDGMENTS 

Very useful discussions with C. Enss, J. Classen, S. 
Hunklinger, S. Ludwig, P. Nalbach, R.O. Pohl, M. The- 
sen, and B. Thimmel are gratefully acknowledged. We 
dedicate this paper to Franz Wegner on the occasion of 
his 60-th birthday, thanking him for numerous discus- 
sions and inspiration throughout the years. 

" El- Mail: kuehn@tphys.uni-heidelberg.de 

E-mail: wuergerQcribxl . u-bordeaux . f r 
^ R.N. Bhatt, in: A. P. Young (cd.) Spin Glasses and Random 

Fields, (World Scientific, Singapore, 1998), p. 225 
^ B.K. Chakrabarti, A. Dutta, and P. Sen, Quantum 

Ising Phases and Transitions in Transverse Ising Models, 

(Springer, Berlin, 1996) 

W. Wu, B. EUman, T.F. Roscnbaum, and G. Acppli, Phys. 
Rev. Lett. 67, 2076 (1991); W. Wu, D. Bitko, T.F. Rosen- 
baum, and G. Aeppli, Phys. Rev. Lett. 71, 1919 (1993) 
* V. Narayanamurti and R.O. Pohl, Rev. Mod. Phys. 42, 201 

(1970) 

^ F. Bridges, Grit. Rev. Solid State Sci. 5, 1 (1975) 



^ A.T. Fiory, Phys. Rev. B 4, 614 (1971); Rev. Sci. Instr. 42, 

930 (1971) 

R. Weis ot al., Phys. Rev. Lett. 75, 2220 (1995); R. Weis 

et al, Ann. Phys. 6, 263 (1997) 
* G. Enss, M. Gaukler, M. Nullmeier, R. Weis, A. Wiirger, 

Phys. Rev. Lett. 78, 370 (1997) 
^ A. Wiirger, R. Weis, M. Gaukler, and C. Enss, Europhys. 

Lett. 33, 533 (1996); C. Enss et al, Phys. Rev. B 53, 12094 

(1996) 

^° A. L. Burin and Yu. Kagan, Physica B 194-196, 393 (1994) 
" A. L. Burin, D. Natelson, D.D. Osheroff, Yu. Kagan, in: 
P. Esquinazi (ed.). Tunneling Systems in Amorphous and 
Crystalline Solids, Springer Berlin New York Heidelberg 
(1998) 

A. Wiirger, Z. f. Phys. B 94, 173 (1994) and 98, 561 (1995) 
A. Wiirger, From Coherent Tunneling to Relaxation, 
Springer Tracts in Modern Physics Vol. 135, Springer 
Berlin Heidelberg New York (1997) 
" N.O. Birge et al., Phys. Rev. B 30, 2306 (1984) 

J.J. Tu and A.J. Sievers, Phys. Rev. Lett. 83, 4077 (1999) 
D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35, 
1792 (1975) 

" R. Kiihn and U. Horstmann, Phys. Rev. Lett. 78, 4067 

(1997) 

1* R. Kiihn and J. Urmann, J. Phys. C 12, 6395 (2000) 
M. Mezard, G. Parisi, and M. A. Virasoro, Spin Glass The- 
ory and Beyond, (World Scientific, Singapore, 1987) 

^° R. Weis, PhD-Thesis, Heidelberg (1996), unpublished 
D. Moy, R.C. Potter, and A.G. Anderson, J. of Low-Temp. 
Phys. 52, 115 (1983) 

W. Kanzig, H.R. Hart, and S. Roberts, Phys. Rev. Lett. 

13, 543 (1964) 
23 L. Viana and A.J. Bray, J. Phys CIS, 3037 (1985) 
2*1. Kanter and H. Sompolinsky, Phys. Rev. B 33, 2073 

(1986) 

P. Cizeau an J.P. Bouchaud, J. Phys. A 26, L187 (1993) 
D.R. Grempel and M.J. Rozenberg, Phys. Rev. Lett. 80, 

389 (1997) 

2^ R.O. Pohl, X. Liu, and R.S. Crandall, Current Opinion in 
SoUd State and Materials Sci. 4, 281 (1999) 

2* U. Horstmann and R. Kiihn, Physica B 263-264, 290 
(1999) 



